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Abstract 

Background: Basic leucine zipper (bZIP) transcription factor gene family is one of the largest and nnost diverse 
families in plants. Current studies have shown that the bZIP proteins regulate numerous growth and developmental 
processes and biotic and abiotic stress responses. Nonetheless, knowledge concerning the specific expression 
patterns and evolutionary history of plant bZIP family members remains very limited. 

Results: We identified 55 bZIP transcription factor-encoding genes in the grapevine {Vitis vinifera) genome, and 
divided them into 10 groups according to the phylogenetic relationship with those in Arabidopsis. The chromosome 
distribution and the collinearity analyses suggest that expansion of the grapevine bZIP (VvbZIP) transcription factor 
family was greatly contributed by the segment/chromosomal duplications, which may be associated with the 
grapevine genome fusion events. Nine intron/exon structural patterns within the bZIP domain and the additional 
conserved motifs were identified among all VvbZIP proteins, and showed a high group-specificity. The predicted 
specificities on DNA-binding domains indicated that some highly conserved amino acid residues exist across each 
major group in the tree of land plant life. The expression patterns of VvbZIP genes across the grapevine gene 
expression atlas, based on microarray technology, suggest that VvbZIP genes are involved in grapevine organ 
development, especially seed development. Expression analysis based on qRT-PCR indicated that VvbZIP genes are 
extensively involved in drought- and heat-responses, with possibly different mechanisms. 

Conclusions: The genome-wide identification, chromosome organization, gene structures, evolutionary and expression 
analyses of grapevine bZIP genes provide an overall insight of this gene family and their potential involvement in 
growth, development and stress responses. This will facilitate further research on the bZIP gene family regarding their 
evolutionary history and biological functions. 

Keywords: bZIP transcription factor family, Grapevine, Gene expression. Drought response, Heat stress response 



Background identified in the plant kingdom [1]. The basic leucine 

Transcription factors (TFs) play important roles in regu- zipper (bZIP) transcription factor family is one of the 

lation and control many crucial biological processes in largest and most diverse families [1,2], The bZIP tran- 

plants. Functional characterization of TFs is essential for scription factors harbor a highly conserved bZIP domain 

understanding transcriptional regulatory networks and composed of two structural features: a basic region 

the biological processes they are involved in. Presently, (N-x7-R/K-x9) for sequence-specific DNA binding, and a 

at least 64 families of transcription factors have been leucine zipper, which is composed of several heptad re- 
peats of Leu or other bulky hydrophobic amino acids, 
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conserved motives, described as R/KxxS/T and S/TxxD/E, 
have been confirmed as sites phosphorylated by the Ca^^ 
independent protein kinase and casein kinase II, respect- 
ively [5]. Proline-rich, glutamine-rich and acidic domains 
are also thought to play roles in transcription activation of 
bZIP genes [6], 

The bZIP gene family has been comprehensively iden- 
tified or predicted in several plant genomes. Seventy-five 
bZIP genes have been found in Arabidopsis {Arabidopsis 
thaliana) [4], 89 in rice {Oryza sativa) [2], 131 in soy- 
bean {Glycine max) [6], 92 in Sorghum [7], and 125 in 
maize {Zea mays) [8]. Among all members in all plant 
species, just a small number of bZIP genes have been 
functionally characterized. Plant bZIP proteins have been 
found to regulate the integration and development of 
many organs and tissues, including seed maturation 
and germination [9,10], floral induction and develop- 
ment [11-18], vascular development [19], embryogen- 
esis [20,21] and photomorphogenesis [22,23]. They 
have also been found to be involved in responses to a 
variety of abiotic/biotic stimuli, such as drought 
[24,25], high salinity [26], heat stress [27], cold stress 
[28,29], pathogen infection [30], and hormone signal- 
ing, such as abscisic acid (ABA) [25,31,32], ethylene 
[33] and light signaling [34]. 

Grapevine {Vitis spp.) is a popular fruit crop that in 
cultivated throughout the world and represents one of 
the most important crops with highly valued products 
such as juice, liquors and wines. The grapevine {Vitis 
vinifera) genome sequence was released in 2007 [35]. 
It was suggested that the grapevine genome was com- 
posed of three ancestral genomes, and it has retained 
more ancient imprint of the evolution. However, the no- 
tion of grapevine as a "palaeo-hexaploid" organism [35,36] 
remains to be confirmed. Since Vitis is sister to all Rosids 
(I and II), studying the evolution of a large gene family in 
grapevine will shed light on our understanding the evolu- 
tion of other plants [36], and will allow the elucidation of 
function and divergence of gene families in grapevine. So 
far, only a limited number of VvbZIPs have been isolated 
and functionally characterized [37]. Among them, two 
VvABFs (homologous of VvbZIPOS and VvbZIP45 in our 
study), isolated from the V vinifera cultivar 'Thompson 
Seedless', were induced by various abiotic stresses including 
drought, salinity, cold and the exogenous hormone ABA 
[38], and VvbZIP45 corresponds to the very recently de- 
scribed ABSCISIC ACID RESPONSE ELEMENT-BINDING 
FACT0R2 {VvABF2\ involved in ABA-dependent grape 
berry ripening [39]. Another gene, VvbZIP23, which was 
isolated from the cultivar 'Mango', was induced by a wide 
spectrum of abiotic stresses, including drought, salt, cold, 
and application of abscisic acid, methyl viologen, sali- 
cylic acid, jasmonic acid, and ethephon [37]. In addition, 
the transcriptional analysis of late ripening stages of 



grapevine berry detected one bZIP transcription factor 
which exhibited lower expression at the harvest stages 
(TH) than that at the 7-days before harvest stage (TH-7) 
[40]. The transcriptomic analysis of grape leaves suggested 
that some VvbZIP transcription factors were involved in 
stress responses [41]. 

In the present study, we report the identification of 55 
bZIP genes from the grapevine genome, and the analysis 
of the gene family with a focus on the evolution and diver- 
gence after multiple duplications in relation to the grape- 
vine genome fusion. We also analyze gene expression of 
VvbZIP genes in different tissues/organs/developmental 
stages, as well as in response to stresses. Our results pro- 
vide a prospective for the evolutionary history and general 
biological involvement of the grapevine bZIP transcription 
factor family. 

Results 

Identification and nomenclature assignment of bZIP 
transcription factors in grapevine 

In order to run a complete search for VvbZIP genes 
in the grapevine genome, all annotated proteins of the lat- 
est version (VI) of the grapevine genome from CRIBI 
(http://genomes.cribi.unipd.it/) were considered. The Hid- 
den Markov Model (HMM) profile of the bZIP domain 
(PF00170) (http://pfam.sanger.ac.uk/) was then employed 
as a query to search the database using the program 
HMMER3.0 with the default E-value. After determin- 
ing the integrity of the bZIP domain using the online 
program SMART (http://smart.embl-heidelberg.de/) with 
an E-value < 0.1 and sequence alignment, 55 non-redundant 
genes were assigned as grapevine bZIP genes, and named 
from VvbZIPOl to VvbZIPSS based on the coordinate 
order on grapevine chromosomes from top to bottom. 
The nomenclature and corresponding information are 
listed in Table 1. 

Phylogenetic analysis and classification of the grapevine 
bZIP gene family 

To investigate the phylogenetic relationship of the bZIP 
gene families in grapevine and Arabidopsis, the amino 
acid sequences of the bZIP domain of the 55 proteins 
from grapevine {Vitis vinifera)^ 72 from Arabidopsis 
{Arabidopsis thaliana) (three genes were no longer sup- 
ported by their annotations comparing to Jakoby et al. 
[4]) and 99 from poplar {Populus trichocarpa) (data from 
Phytozome v9.1) were used to construct a phylogenetic 
tree. Three methods, Neighbor-Joining (NJ), Minimal Evo- 
lution (ME) and Maximum Parsimony (MP), yielded nearly 
identical phylogenetic trees (data not shown), therefore, 
only the NJ tree was used for further analyses (Figure 1). 
Based on the phylogenetic tree, the VvbZIPs can generally 
be divided into the same 10 subgroups as those in Arabi- 
dopsis. However, AtbZIP62 1 VvbZIP26/PotrL006G 114600.1, 



Table 1 The identified VvbZIP genes and their related information 



Gene 
name 


Gene ID 


Group 


Chr. 


Locus 


Chain Protein 
Length 


Amount 
of intron 


Intron pattern 
within bZIP domian 


AtbZIP Ortholog (Syntenic orthologs were denoted by 
overstriking and italics) 


VvbZIPOl 


VIT_01s0011g03230 


J 


1 


2941266 


2946275 


+ 453 


10 


b 


46, 20, 45, 26, 65, 21, 57, 47, 22, 50 


VvbZIP02 


VIT_01s0010g00930 


C 


1 


16444033 

1 \J III \J ^ ^ 


1 6444682 


195 


0 


h 


42, 58, 48, 43, 3, 8, 70, 2, 44, 53, 7, 1 1, 41, 5, 68, 6, 16, 4 


VvbZIPOS 


VIT_02s0025g01020 


D 


2 


989255 


996186 


+ 398 


10 


0 


55(GBF3), 54(GBF2), 16, 68, 41 


VvbZIP04 


VIT_02s0012g02250 


D 


2 


9137736 


9153313 


413 


11 


a 


55, 54, 76, 68, 41, 44, 42, 43, 58, 3, 70, 8, 48 


VvbZIPOS 


VIT_03s0038g00860 


F 


3 


689247 


693308 


+ 676 


5 


c 


29,30, 18, 52, 69, 51, 59, 31, 33, 61, 34, 74 


VvbZIP06 


VIT_03s0038g02420 


1 


3 


1679507 


1683223 


+ 271 


1 


h 


23, 19, 24 


VvbZlPO/ 


VIT_03s0038g04450 


C 


3 


3227416 


3228660 


157 


0 


h 


44, 1 1, 2, 53, 48, 3, 42, 43, 58, 8, 1, 6, 7, 5, 68, 70 


VvbZIPOS 


VIT_03s0063g00310 


A 


3 


^88707^ 


3890514 


409 


4 


0 


37(ABF3), 36, 38(ABF4/AREB2), 35(ABF1), 15, 66, 12, 67, 39, 14, 13, 40 


VvbZIP09 


VIT_04s0008g02750 


B 


4 


2284121 


22871 15 


349 


5 


d 


9, 63, 1 0, 2, 53, 42, 25, 44, 58, 3, 48, 68, 1 1 


VvbZIP 10 


VIT_04s0008g05210 


G 


4 


4716191 


471 9252 


+ 169 


3 


f 


56(HY5), 64 


VvbZIP 11 


VIT_04s0069g01150 


A 


4 


9841491 


9842587 


+ 57 


1 


a 


66, 12, 37, 36, 35, 38, 39 


VvbZIP 12 


VIT_04s0023g01360 


D 


4 


1 7832626 


178639S2 


+ 364 


10 


0 


47, 16, 68, 55, 54, 42, 70, 3, 58, 43, 48, 8, 44, 18 


VvbZIP 13 


VIT_04s0023g02430 


C 


4 


18983616 


1 8984767 


+ 249 


0 


h 


44, 1 1, 2, 53, 48, 42, 58, 3, 43, 8, 7, 6, 1, 5 


VvbZIPU 


VIT_05s0077g01140 


C 


5 


878161 


879142 


154 


0 


h 


53, 2, 44, 11,58, 3, 5, 48, 7, 1 


VvbZIPlS 


VIT_05s0020g01090 


G 


5 


7860048 


2863747 


+ 210 


3 


h 


64(HYH), 56 


VvbZIP 16 


VIT_05s0102g01120 


F 


5 




23231 749 


425 


3 


c 


69, 59, 18, 52, 51, 29, 30, 61, 33, 31, 34, 74 


VvbZIPU 


VIT_06s0004g08070 


F 


6 


883S674 


8841 OSl 


342 


3 


c 


78,52, 59, 69, 51, 29, 30, 33, 61, 31, 74, 34 


VvbZIPlS 


VIT_06s0009g01790 


A 


6 


1 3809459 


1 381 S62S 


248 


2 


0 


66, 12, 35, 38, 36, 67, 39, 40, 14, 13, 37, 15 


VvbZIP 19 


VIT_06s0080g00340 


A 


6 


20262007 


20263157 


+ 324 


2 


a 


67(DPBF2), 39, 66, 12, 35, 36, 38, 15, 37, 40 


VvbZIP20 


VIT_06s0080g00360 


J 


6 


20304714 


2031 1349 


500 


10 


g 


27, 65, 45, 20, 26, 46, 22, 50, 47, 57 


VvbZIP21 


VIT_07s0141g00170 


B 


7 


1 16657 


12681 1 


452 


5 


d 


25, 10, 63,9 


VvbZIP22 


VIT_07s0005g01450 


C 


7 


3987031 


39881 20 


+ 145 


0 


h 


53, 44, 2, 1 1 , 42, 3, 6, 48, 7, 58, 43, 8, 1 6, 1 , 68 


VvbZIP23 


VIT_07s0031g01320 


J 


7 


1 7404970 


17410595 


+ 349 


7 


b 


47, 57, 50. 22, 20, 45, 26, 21, 46, 65 


VvbZIP24 


VIT_08s0040g00870 


H 


8 


1 1857207 


1 1860728 


819 


1 


h 


1 7, 28, 49 


VvbZIP25 


VIT_08s0007g03420 


A 


8 


1 73721 00 

1 / ^ / Z_ 1 \J\J 


1 73780S7 


400 


0 


a 


39(ABI5), 35, 67, 66, 38, 37, 36, 12, 13, 15, 40, 14 


VvbZIP26 


VIT_08s0007g03640 


UC 


8 


\ / j/yoD4 


1 ~/^Q(^ 1 QQ 

1 / joD 1 by 


+ 496 


5 


0 


62 


VvbZIP27 


VIT_08s0007g05170 


J 


8 


19088004 


19098613 


451 


10 


b 


45, 20, 26, 46, 21, 65, 57, 50, 22, 47 


VvbZIP2S 


VIT_08s0007g06160 


J 


8 


19963578 


19968190 


491 


10 


b 


65, 21, 45, 20, 26, 46, 22, 57, 50, 47 


VvbZIP29 


VIT_12s0028g02590 


E 


12 


3366293 


3367561 


295 


3 


e 


61, 34, 51, 69, 18, 52, 59, 30, 29, 72, 33, 31 


VvbZIPSO 


VIT_12s0055g00420 


A 


12 


13276799 


13300251 


+ 299 


3 


0 


13, 40, 12, 66, 67, 35, 37, 38, 36, 39, 14 


VvbZIPSl 


VIT_12s0034g00110 


A 


12 


15485471 


15491224 


+ 99 


3 


a 


66, 12, 37, 36, 35, 38, 39 


VvbZIP32 


VIT_12s0035g00620 


E 


12 


20038502 


20045379 


+ 374 


3 


e 


61,34, 69, 59, 18,51,29, 52, 30 



3 

Q. O 



3 

^ 00 



O 
-h 

00 
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VvbZIP33 


VIT_ 


.13s0067g02900 


F 


13 


1562503 


1569804 


+ 


359 


3 


VvbZIP34 


VIT_ 


.13s0175g00120 


A 


13 


15049338 


15072969 


+ 


325 


2 


VvbZIP35 


VIT_ 


.13s0084g00660 


J 


13 


19780668 


19805654 


- 


469 


10 


VvbZIP36 


VIT_ 


.13s0158g00380 


1 


13 


21283084 


21284766 


- 


286 


1 


VvbZIP37 


VIT_ 


.14s0060g01210 


C 


14 


959578 


961198 


- 


145 


0 


VvbZIP38 


VIT_ 


.14s0030g02200 


B 


14 


7402807 


7408597 


+ 


423 


5 


VvbZIP39 


VIT_ 


.14s0083g00700 


C 


14 


22839619 


22840360 


- 


198 


0 


VvbZIP40 


VIT_ 


.15s0046g01440 


D 


15 


18416554 


18423462 


- 


430 


11 


VvbZIP41 


VIT_ 


.18s0122g00500 


UC 


18 


420616 


423564 


+ 


322 


2 


VvbZIP42 


VIT_ 


.18s0001g04470 


J 


18 


3861987 


3867762 


- 


363 


7 


VvbZIP43 


VIT_ 


.18s0001g04500 


A 


18 


3875315 


3878509 


+ 


311 


2 


VvbZIP44 


VIT_ 


.18s0001g08710 


C 


18 


7185268 


7185664 


- 


81 


0 


VvbZIP45 


VIT_ 


.18s0001g 10450 


A 


18 


8761978 


8769742 


- 


447 


3 


VvbZIP46 


VIT_ 


.18s0001gl2120 


D 


18 


10309606 


10313842 


+ 


396 


10 


VvbZIP47 


VIT_ 


.18s0001g 13040 


C 


18 


11138196 


11139020 


+ 


154 


0 


VvbZIP48 


VIT_ 


.18s0001g 13740 


F 


18 


1 1 745490 


11746413 


- 


187 


3 


VvbZIP49 


VIT_ 


.18s0001g 14890 


A 


18 


12936974 


12937903 


- 


214 


2 


VvbZIPSO 


VIT_ 


.18s0076g00330 


F 


18 


16209694 


16218061 


+ 


350 


3 


VvbZIP 51 


VIT_ 


.18s0072g00470 


A 


18 


19680318 


19687522 


+ 


275 


3 


VvbZIP52 


VIT_ 


.19s0014g01780 


E 


19 


1945207 


1947498 




309 


3 


VvbZIP53 


VIT_ 


.19s0015g01020 


A 


19 


9101591 


9112651 




251 


2 


VvbZIP54 


VIT_ 


.18s0001g03010 


E 


18_random 


3068050 


3070292 




272 


3 


VvbZIP55 


VIT_ 


00s0541g00020 


C 


Un 


31593865 


31594595 




194 


0 



Note: AtbZIP Orthologs were identified by local BLAST search with E-value < 1E-10 and the identification of 



c 18, 52, 59, 69, 51, 29, 30, 33, 61, 31, 34, 74, 72 

o 66, 12,35,67,38,37,36,39, 15, 14 

b 45, 26, 20, 46, 21, 65, 57, 22, 50, 47 

h 23,19,24 

h 53, 1 1 , 2, 44, 1 , 48, 58, 3, 8, 7, 42, 43, 70 

d 63.10.25.9 

h 42,43, 58, 3,48, 8, 70, 44, 2, 53, 6, 1 1 , 41 , 7, 5, 68 

o 55(GBF3), 54(GBF2), 1 6, 68, 41 , 3, 42, 43, 70, 58 

f 60 

b 57, 47, 50, 22, 26, 20, 65, 45, 21 , 46 

o 1 2, 66, 35, 36, 39, 38, 37, 1 4, 1 5, 67, 40, 1 3 

h 7, 6, 4, 5, 2, 3 

o 36(ABF2/AREB1), 35, 37(ABF3), 38, 66, 12, 39, 75, 67 

o 41 , 1 6, 68, 55, 54, 43, 58, 42, 3, 48, 70, 8 

h 44, 1 1 , 2, 53, 48, 3, 42, 58, 43, 1 , 7, 6, 8, 5 

c 30, 29, 51, 69, 52, 18, 59, 33, 31, 61, 34 

o 14(FD), 66, 1 2, 37, 36, 39, 38, 35, 27 

c 51, 18, 69, 59, 52, 29, 30, 31, 61, 33, 74, 34 

o 13,40,66, 12,39,36 

e 61, 34, 51, 18, 52, 69, 59, 30, 29, 72, 33, 31 

g 36,12 

e 61,34,69,51,59 

h 6, 5, 7, 44, 42, 58, 2, 3, 1 1 , 48, 53, 43, 70, 8, 1 , 4 



lie orthologs (highlight with italics) was performed by using the software MCScanX. 
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„ . AlbZIPI 
^ Potn.OOSG 106700 1 
? Potri.010G142900.1 
AtbZIPS 

-5 AtbZIP4 
C« AtbZIP7 
^ MbZIPe 



A (\IAra) 

Figure 1 The phylogenetic tree of grapevine and Arabidopsis bZIP genes. The 55 grapevine, 72 Arabidopsis and 99 poplar bZIP 
domain protein sequences were aligned by Clustal X 1.83 and the phylogenetic tree was constructed using MGA5.0 by the 
Neighbor-Joining (NJ) method. The Bootstrap value was 1,000 replicates. The gene names of the characterized bZIP genes are 
shown in parenthesis after the gene numbers, and which were referred to previous research [2,4]. A{A/Aro), B(C//\ra) and so on is presented 
as the subfamily names of bZIP gene family in grapevine (capital letter outside the parenthesis) and Arabidopsis (capital letter in parenthesis 
and Ara represent for Arabidopsis), UC, unclassed genes or clades. The colored branch indicates the different subfamily and the black 
represent for UC. 



AtbZIP60/VvbZIP41/PotrL00SG257900.h and AtbZIP72/ 
PotrL009G075000J (UC, Figure 1) form three small unique 
clades in the phylogenetic tree and may have had independ- 
ent evolutionary trajectories from others clades. In addition, 
AtbZIP74, AtbZIPSS and AtbZIPSl and PotrL006G237500J/ 
PotrL018G04S100J/PotrL001G077900J/PotrL001G076S00J/ 
PotrL003G220000J/PotrW01G004900.h were distinguished 
from other bZIP genes as they formed individual clades 
without any other bZIP genes, suggesting that these indi- 
vidual clades may be specific to Arabidopsis and poplar, 
respectively. This result, different from the former classifi- 
cation in Arabidopsis [4], was supported by Nijhawan 
et al. and Wei et al. [2,8] and sustained by the analyses of 
the gene structures, conserved motifs and DNA-binding 
site specificity of VvbZIP transcription factors in the 
present work. 



Chromosomal location and expansion patterns of VvbZIP 
genes 

The Vitaceae evolved about 60 million years ago (MYA) 
and is the earliest diverging lineage of rosids [36,42]. The 
extant grapevine genome was considered to be formed by 
fusion of three ancestral genomes following post-fusion 
shuffling [36] which has not undergone recent genome 
duplication following the y WGD. The "palaeo-hexaploid" 
nature of the genome is supported by many triplicate 
blocks in the genome, as indicated by the same colors of 
homeologous blocks adapted from Jaillon et al. [35], and it 
was hypothesized that the ancestral hexaploid genome of 
V. vinifera was derived from an allopolyploidization of 
two ancestral genomes, which were named as Va and Vc, 
respectively [36]. Chromosome triplets consist of two Va 
and one Vc chromosomes, as expected from the tetraploid 



Liu et at. BMC Genomics 2014, 15:281 
httpy/www.biomedcentral.com/l 471 -21 64/1 5/281 



Page 6 of 18 



chr1 
23.0M 



chr2 
18.7M 

_VvbZIP03 



chr4 
23.8M 



_VvbZIP01 



chr5 chr6 
25.0M 21. 5M 

.^vbZIP14 
_wbzTF 








l_VvbZIP41 


_VvbZIP52 






^VvbZIP42 
\VvbZIP43 








_VvbZIP44 

_VvbZIP45 
_VvbZIP46 


_VvbZIP53 






>\VvbZIP47 








VVvbZIP48 
\VvbZIP49 

\VvbZIP50 






VvbZIP51 





Figure 2 Distributions of VvbZIP genes in grapevine Chromosomes. Fifty-three VvbZIP genes were mapped to the 14 linkage groups. 
VvbZIP54 was located on 18 random chromosomes and VvbZIP55 was assigned to the Un chromosome. Five chromosomes (chr9, 10, 1 1, 16 and 17) 
have no grapevine bZIP gene, The schematic representation of paralogous regions was adapted from Jaillon et al., with the same color [35]. Acording 
to the former hypothesis [36], the ancestral hexaploid genome of V. vinifero was derived from an allopolyploidization of two ancestral genomes: Va 
and Vc, which are denoted on the right of the figure, and Un is unassigned chromosomes. The chr names in red underline the chromosomes that 
harbored no bZIP genes. The gene names with collinearity were performed with the some color of chromosome paralogous region. 



and diploid conditions of the two component genomes 
[36]. To search for imprints of the grapevine genome evo- 
lution from the perspective of gene families, the location 
of each VvbZIP gene was given a diagrammatic represen- 
tation based on the annotation of the 12x VI assembly 
grapevine genome available at CRIBI (Figure 2). The re- 
sults showed that the distribution of VvbZIP genes is 
not even. VvbZIP54 and VvbZIPSS are located on the 18 
Random and Un chromosome, respectively. Among the 
other 53 VvbZIP genes, 11 (20.75%) genes were situated 
on chromosome 18; 5 (9.43%) VvbZIP genes were each lo- 
cated in the chromosome 4 and 8; 4 VvbZIP genes posi- 
tioned each in chromosomes 3, 6, 12 and 13; 3 VvbZIP 
genes located each in chromosomes 5, 7 and 14; 2 VvbZIP 
genes each in chromosomes 1, 2 and 19; only one gene in 
chromosome 15; and no VvbZIP genes are located in 
chromosomes 9, 10, 11, 16 and 17. 

Furthermore, we analyzed the tandem arrangement 
and the collinearity by using MCScanX program and 
mapped the genes to chromosomal locations (Figure 2), 
and found 14 homologous blocks containing bZIP genes 
with collinearity that may have resulted from genome fu- 
sion or WGD and none of the VvbZIP genes was found 



to be arranged in tandem. Interestingly, all of those 
collinearity events occurred within same homologous 
blocks, and only two triplicate bZIP genes, VvbZIP?/ 
13/47 and VvbZIP 14/22/37, were identified, indicating 
that the expansion patterns of bZIP gene family par- 
tially support the triplicate state of the grapevine gen- 
ome and combined with further genome fusion or WGD 
events. Further, the rest of the collinearity events were just 
found within two of the homologous triplets, and the third 
homologous triplet, such as chromosome 9, 10, 11, 16 
and 17, mostly harbored no VvbZIP genes, which 
seem to belong to Vc genome of grapevine described 
by Malacarne et al. [36], supposing that VvbZIP genes 
in the heterogenous Vc genome may have undergone 
extensive loss and shuffling after genome fusion. 

Gene structure of VvbZIP genes 

Since the intron/exon organizations and intron types 
and numbers are typical imprints of the evolution within 
some gene families [2,43,44], the VvbZIP gene structures 
were examined to obtain further insight into their evolu- 
tionary trajectory. The results show that exon/intron struc- 
tures of VvbZIPs were highly conserved within subfamilies. 
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Genes that clustered together generally possessed a similar 
gene structure, especially the exon number and length 
(Additional file la). Interestingly, no intron containing 
genes were found in group C, which accounted for 18.2% 
(10/55) of the total VvbZIP genes. Furthermore, the num- 
ber of exons in each group (Additional file lb) was uneven. 
The numbers of exons in group B, C, E, G and I were 
much more conserved than that in those in the group A 
(varying from one to 5). Furthermore, the positions and 
phases of introns in the bZIP domain (composing of basic 
region, hinge and leucine zipper) were analyzed (Figure 3). 
Nine intron/exon structural patterns in the VvbZIP gene 
family were identified based on the intron position, 
number and splicing phase (details were represented in 
Additional file 2), similar to those in rice and maize [2,8]. 
We named these patterns as a, b, c, d, e, f, g,h and /. The 
patterns a, b and / in VvbZIP correspond to the patterns 
ay b and /in maize and the patterns b, c and g in rice, re- 
spectively. The patterns c, d, e and /in VvbZIP jointly cor- 
respond to the pattern c in maize or the pattern a in rice, 
because they have the same position of intron in the basic 
and hinge regions. 

In addition, each subfamily of VvbZIPs mostly shared 
the same intron/exon structural pattern. As illustrated in 
Figure 3, Pattern a has just one intron in phase 0 (PO) in 
the leucine zipper which inserts between the amino 
acids Gin and Ala, includes members of group A with- 
out VvbZIP20 and all genes of the group D. Pattern b 
has two introns in phase 0 (PO) in the basic and hinge 
regions, respectively, and includes the genes of group J 
without VvbZIP20, The pattern c, d, e and / share a 
common intron insert with phase 2 (P2) in the basic re- 
gion, but for each pattern, differences in intron inserts 
exist in the leucine zipper region. Pattern c, d and e have 



another intron insert in phase 0 (PO) in the leucine zip- 
per region at a different position, whereas the pattern / 
has no intron insert in the leucine zipper region thus 
different from c, d and e; and the pattern c, d and e are 
completely corresponded to the group F, B, E, respect- 
ively. Pattern g has only one intron in phase 1 (PI) in 
the basic region and pattern h has one intron in phase 0 
(PO) in the leucine zipper region; both patterns g and h 
contained just one member of VvbZIP gene each from 
group J and group A, respectively. Pattern / lacks any in- 
tron in the basic region and leucine zipper region, and 
this pattern contains all genes of group C (10 genes), H 
(one gene), group I (two genes) and VvbZIP 15 in the 
group G. In summary, the intron position and splicing 
phase in the basic region and the leucine zipper region 
of grapevine bZIP gene family have been highly con- 
served during the course of evolution in each group, 
which is concordant with the early reports [2,8] . 

Furthermore, considering the patterns of bZIP gene 
family of rice, maize and grape, we find that although 
there are some differences of the splicing phase and the 
interrupted amino acids, the positions of the intron in- 
sert are highly conserved. This phenomenon suggested 
that the position of the intron insert in the conserved 
domain is important for the evolution and divergence of 
function in the gene family. 

Conserved structural features of VvbZIP transcription 
factors 

The bZIP domain is the core of the bZIP gene family 
members, which preferentially binds to the promoter of 
their downstream target genes on a specific cis -element 
(e.g. ABREs). However, the functional diversity of the 
bZIP gene family may also be contributed by additional 



Pattern of 
introns 




RfRQRR|VllKNRESAARSRARKCjAYTMELEA^VAKLI^EKNEEL 

N— X(7) R/K X(9) L — X(6) — L X(6) — L 

p 



4i- 



Numbers of 
VvbZIP 
genes 

18 

6 

» 6 
3 
4 
2 
1 
1 

14 



VvbZIP genes 
in each pattern 

Group A without 
VvbZIP53, Group D 

Group J without VvbZIP20 

Group F 

Group B 

Group E 

VvbZIP 10(G), VvbZIP41(UC) 

VvbZIP20(J) 

VvbZIP53(A) 

Group C, Group H, Group I and 
VvbZIP15(G) 



Figure 3 A schematic diagram of Intron patterns In the grapevine bZIP gene family. The primary structure of the bZIP domain is shown at 
the top of the picture. The basic region, hinge and leucine zipper were shaded in crimson, blue and green, respectively. The highly conserved 
residues are highlighted with red and light green boxes. An example of the consensus sequence is given below, and the corresponding positions 
are shown at the top of the bar [45]. Leucines are sometimes replaced by He, Val, Phe or Met. The black vertical lines in the bar indicated the 
intron position within the bZIP domain and PO, PI and P2 stand for the intron phase 0, phase 1 and phase 2, respectively. The details of intron 
position within bZIP domain of the VvbZIP genes were shown at Additional file 2. 
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domains in the bZIP proteins [46] . In total, 20 conserved 
motifs, including the bZIP domain, were identified in the 
VvbZIP proteins and their multilevel consensus amino 
acid sequences of motifs are listed in Additional file 3. 
The motif distribution corresponding to the phylogenetic 
tree of VvbZIP gene family is shown in Figure 4a,b. Motif 
1 is the basic region and the hinge of the bZIP domain; 
whilst motif 5 (typically), 3 and 9 represent the variable 
motifs in the leucine zipper region across the bZIP gene 
family. For example, motif 3 only appeared in the group E 
and F; motif 9 appeared in the group J exclusively; and 
group C and D mostly harbored a double motif 5. 
Further, the number of Leucine units ranged from 3 to 
11, and some of them were interrupted by one or two 
other units which are conserved with the same number of 
amino acids. Finally, nine gradients (I to IX) were defined 
(Additional file 4), and some of them showed group speci- 
ficity, such as the gradient I which harbored the most 
members of group A, and gradient IX which harbored the 



most members of group C. Some typical cases and their 
corresponding motifs are shown in Figure 4c and may 
provide a greater understanding of the divergence and the 
roles of leucine zipper region in the bZIP gene family. 

Beside the bZIP domain, each group of VvbZIP shares 
some motifs. For example, most members of the group 
A possess motifs 6, 7, 8 and 20; group B and C share 
motif 11; all members of the group D specifically harbor 
the motifs 10 and 13; group E shares motif 19 with 
group F which also possesses motif 18 exclusively. Motif 
16 appeared in the group I, motifs 2, 4, 9, 14 and 15 are 
present in group J exclusively, whereas group G and H 
lack any additional motif. 

Previous evidence suggested that most members of 
the group A of the bZIP gene family are involved in ABA- 
responsive signal-transduction pathway and abiotic stress 
response [4,25,47,48]. Typically, ABF/AREBs are master 
transcription factors that regulate ABRE-dependent ABA 
signaling involved in drought stress tolerance, but their 
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Figure 4 Conserved motifs of VvbZIP proteins. All conserved motifs of the VvbZIP proteins were identified by the MEME program. The protein 
structures (b) of VvbZIPs based on the presence of conserved motifs were arranged corresponding to the phylogenetic tree (a) of the VvbZIP 
proteins. Different motifs are highlighted with different colored boxes with numbers 1 to 20. Among them, the motif 1, 5, 9 and 13 mainly 
represented variations in the bZIP domain and the unusual motifs is shown in white. The sequence details of these conserved motifs are given in 
Additional file 3. The 3D structure and the variation of Leucine zipper region of VvbZIP proteins in length were described in (c). 
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activation requires posttranscriptional modification with 
ABA-dependent phosphorylation [25]. In the present 
study, two conserved sites, R/KxxS/T and S/TxxD/E, were 
verified as the key phosphorylation sites in VvbZIP pro- 
teins which are phosphorylated by a Ca^^ independent 
protein kinase and casein kinase II, respectively [5,32]. 
Furthermore, ABA-dependent multiple phosphorylation 
sites are also required for their activation. In our study 
(Figure 4b), motif 6, 7, 8 and 20 appeared in the group A 
VvbZIP exclusively and they widely represent the potential 
phosphorylation sites described above. Motif 6 harbors 
the R/KxxS/T and motif 7 harbors S/TxxD/E, presented 
as RQ[GNS]S, T[LF]DE and R[QE]PT, TLE[DE], re- 
spectively. Motif 8 and 20 contain just one of these sites, 
presented as TV[DE]E and R[TV][SL][ST], respectively. 
Moreover, these potential phosphorylation sites were also 
found in other groups of the VvbZIP gene family, pre- 
sented as RQQT/ RALS in motif 4, [RK][LV]Q[SM] in 
motif 11, T[DI][VAM][DEQ] in motif 15, and R[SA]SS in 
motif 19, respectively. In addition, other conserved motifs 
in bZIPs that may function as transcriptional activators, 
have also been found, including proline-rich, glutamine- 
rich, and acidic domains. 

Predicted DNA-binding site specificity of VvbZIP 
transcription factors 

As one of the biggest transcription factor families in 
eukaryotes, bZIPs owe their DNA-binding ability and 
specificity to the bZIP domain. We aligned the amino 
acid sequences of each subfamily and compared them 
with those in Arabidopsis, rice and other important 
species in the plant kingdom, including Physcomitrella 
patens (moss), Selaginella moellendorffii, Amborella 
trichopoda (Additional file 5). The results show that 
some amino acid residues are highly conserved within 
each group across species in the tree of land plant life, in- 
dicating that those amino acid residues play crucial roles 
in the DNA-binding ability of the bZIP domain and have 
undergone purifying selection during their evolution. 
Experiments have shown that asparagine (numbered 
as -18), arginine (numbered as -10) and leucines (num- 
bered as +1, +8, +15), were key and invariant sites of 
the bZIP domain, and their replacement causes new 
binding specificities [8,45,49]. In group F and in VvbZIP26 
which seemingly belonged to group D in the phylo- 
genetic tree (Figure 4a), the arginine site was replaced 
by Lys (just occurred in group F) and He, respectively, 
and is consistent with its subfamily members in other 
species (UC, Additional file 5). These replacements 
may have contributed to the emergence of a new 
binding specificity and independent course of evolu- 
tion of this subfamily. In addition, some leucine sites 
were also replaced by other amino acids (for example, 
the third leucine site (+15) of most members in 



group C were replaced by He, and the second leucine 
site (+8) were mostly replaced by Met) which may 
have led to change in the binding specificity of these 
bZIP domains as well. 

Expression analysis of VvbZIP transcription factors at 
different developmental stages in specific organs in 
grapevine 

Evidence has suggested that bZIP transcription factors 
are widely involved in the integration and development 
of many organs and tissues, such as seed maturation and 
germination [9], floral induction and development [11-18]. 
To obtain more insight into the temporal and spatial tran- 
scription patterns of grapevine bZIP gene family during 
grapevine development, a heatmap of the global tran- 
scription profile of the VvbZIP family was performed 
based on microarray data of the recent grape gene ex- 
pression atlas on 54 samples (the sample information is 
described in Additional file 6) [50]. Fluorescence in- 
tensity values were log2-based to narrow the dynamic 
range of fluorescence intensities, thus giving a better 
appreciation of small differences in gene expression. All 
members of the grapevine bZIP gene family were sur- 
veyed. As illustrated in Figure 5, the VvbZIP gene family 
has a broad expression pattern across various develop- 
mental stages of grape organs and tissues. The VvbZIP4S, 
VvbZIP41, VvbZIP24 and VvbZIPOS were abundantly 
expressed in all examined organs and tissues. In con- 
trast, VvbZIP39, 43, 55, SI and 44 were expressed at 
relatively low level. Moreover, VvbZIPSS, and VvbZIP 1 9 
were mainly expressed in seed at the stage of post-fruit 
set while VvbZIP25 expression increased during the 
whole seed development. Notably, those three genes 
belonged to a subclade of the group A in the phylogenetic 
tree (Figure 1) harboring Arabidopsis homologous genes 
{AtbZIP67/DPBF2, AtbZIP39/ABIS and AtbZIPlS), which 
have been confirmed to play an important role in ABA- 
mediated seed development, germination, and embryo 
maturation [51,52]. This suggests that VvbZIP 19, 25 and 
53 may have functions similar to their orthologous 
Arabidopsis genes during seed development in grape- 
vine. Interestingly, the group E of VvbZIP showed pe- 
culiar expression in flower organs. For example, the 
VvbZIP29 and 32 were mainly expressed in pollen 
and the VvbZIP54 was highly expressed in Flower-FB, 
Flower-F, stamen and pollen, and the last member of 
group E, VvbZIP52y was highly expressed in the carpel but 
was relatively low in pollen, among others, indicating that 
the whole group E of grapevine bZIP gene family may be 
involved in flower development. 

To survey the co-expression relationship between VvbZIP 
genes, Pearsons correlation coefficients (PCC value) were 
calculated among bZIP genes based on the microarray 
data of 54 samples (Additional file 7), and showed that 
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Figure 5 Expression profiles of VvbZIP genes in 54 samples including green and woody tissues and organs at different developmental 
stages and some specialized tissues. Fluorescence intensities values were log2-based. Genes were listed according to those in the phylogenetic 
tree and samples names are based on the gene expression atlas of grapevine [50]. 



significant correlations existed between 213 pairs of 
VvbZIP genes at the 0.01 level (2-tailed), including 
154 (72.3%) pairs showed significant positive correlation 
(PCC value ^0.35) and 59 (27.7%) pairs showed sig- 
nificant negative correlation (PCC value ^ - 0.35). Among 
them, just four pairs of genes, VvbZIP 19 (Group A)/53 
(Group A), VvbZIP2S (Group A)/28 (Group J), VvbZIPll 
(Group K)I31 (Group A) and VvbZIP29 (Group E)/54 
(Group E), showed strong positive correlations with PCC 
value of > 0.9, and seemed to belong to the same group, in- 
dicating that they probably contribute to some redundant 
function or are involved in the same regulation network of 
biological process. 

In addition, we carefully examined the differentially 
expressed VvbZIP genes in the development of various 
organs (Additional file 8), suggesting that the VvbZIP 
genes are involved in seed development (e.g., VvbZIP25, 
14, 20, 53 19, 28, 7, 47 and 52), leaf and stem develop- 
ment (e.g., VvbZIP20, 23, 44, 14 and 17), berry matur- 
ity (e.g., VvbZIP2, 44, 42 and 49) and bud dormancy. 
Interestingly, VvbZIP 14, 44, 20 and 47 were involved 
in the development of more than one organ, indicating 



that they may participate in some basic biological 
processes. 

Analysis of VvbZIP gene expressions by quantitative real 
time PGR 

Quantitative real-time PCR (qRT-PCR) was performed 
on four different organs (stem, root, leaf and shoot tip) 
of Pinot Noir PN40024 grapevine plants for VvbZIP 
genes, and the results obtained were compared with 
microarray data of the cultivar 'Corvina [50] with the 
same organs at the corresponding phonological phases 
(green stem (Stem-G), root (Root), mature leaf (Leaf-FS) 
and bud after-burst (Bud-AB)) (Figure 6). The results 
show that the expression profiles of VvbZIP genes from 
qRT-PCR were mostly in agreement with the microarray 
data based on calculation of PCC value. This means that 
similar spatiotemporal expression of members of the 
VvbZIP gene family is generally conserved among different 
cultivars like Corvina and Pinot Noir, grown in very differ- 
ent conditions, i.e. in the field and in vitro, respectively. 
Our results provide evidence for selecting candidate genes 
for further characterization in tissue/organ specificity. 
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Figure 6 The comparison between quantitative real-time RT-PCR data and Microarray data. The relative expressions of 24 selected VvbZIP 
genes were performed with quantitative real-time RT-PCR (q RT-PCR). The relative expression of stem/green stem was set up as 1. The y-axis 
indicates the folds of gene expression relative to stem. The expression in significant agreement between qRT-PCR data and Microarray 
data were marked with stars based on PCC value in 0.01 (two stars). Characters on the x-axis indicate various tissues, S, stem/green stem; R, root/root; 
L, leaf/ mature leaf; ST, shoot tip/ bud after-burst. The black bars denote the qRT-PCR data, and the gray bar denote the microarray data. 



Expression of VvbZIP genes under drought and heat 
stress 

Plants and crops are frequently challenged by abiotic 
stress, such as drought stress and high temperature. Some 
evidence has suggested that bZIP proteins are widely 
involved in signaling and responses to abiotic/biotic stim- 
uli [2,26,30,53], but limited information was available on 
bZIP involvement on stress responses in grapevine. In the 
present study, long-term drought condition (0, 4, 8 and 
12-day) and short-term heat stress (0, 2, 5 hours and 
24-hour recovery) on grapevine seedling treatment had 
marked effect on the expression profile of 47 out of 55 
(85.5%) VvbZIP genes for which we were able to design 
primers using for qRT-PCR analysis (Figure 7). As illus- 
trated in Figure 7, both heatmaps (responses to drought 
and heat) could be divided into four clusters. For the 
drought treatment, cluster 1 contains 6 (13.3%) members 
of detectable VvbZIP genes, which were widely down- 
regulated by drought treatment. Cluster 2 contains genes 



highly induced at 4 days after drought treatment, but 
down-regulated at 8 days and 12 days after drought treat- 
ment. Cluster 3 mainly consists of genes with the highest 
expression level at 8 days after drought treatment. 
Cluster 4 has 15 (33.3%) members which were mainly 
up-regulated with the increased degree of the drought 
treatment, and particularly, VvbZIP44, 45, 14, 55 and 
03 were highly induced at 12 days after drought treat- 
ment. Overall, 25 VvbZIP genes were up-regulated at 
least two fold after drought treatment relative to the 
control and ranged from 2 to 25 fold, 7 VvbZIP genes 
were down-regulated at least two fold after drought 
treatment relative to the control and ranged from 2 
to 239 fold. 

For the heat stress treatment and recovery, cluster 1 has 
three (6.7%) members of 45 detectable VvbZIP genes, 
which were widely down-regulated after 45°C treatment 
and recovery. The members of cluster 2 were mainly 
down- regulated after 24 hour of recovery. The members 
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Figure 7 Quantitative RT-PCR analysis of the VvbZIP gene expression in grapevine leaves (sequenced genotype, PN40024) in response 
to drought stress and 45°C heat stress. The normalized relative expression level of the 0 d treatment timepoint was set up as 1. The y-axis 
indicates the relative folds of gene expression compared with the control. Leaves collected at 0, 4, 8 and 12 d post-drought or at 0, 2, 5 h post 
heat stress and 24 h recovery. Means are calculated from technical triplicate qRT-PCR measurements within three biological replicates. 



of cluster 3 were mainly down-regulated after 45°C treat- 
ment for 2 or 5 hours, but were up-regulated after recov- 
ery for 24 hours. Cluster 4 expressed highest across all 
points of heat stress treatment. Overall, 14 VvbZIP genes 
were up-regulated at least two folds after 45°C treatment 
or recovery, the fold change ranged from 2 to 7; 7 VvbZIP 
genes were down-regulated at least two fold and ranged 
from 2 to 7. 

Interestingly, some up-regulated genes typically in- 
duced by drought were down-regulated during the heat 
stress or in the recovery process, such as VvbZIPOl, 39, 
22 and 55. On the contrary, some genes down-regulated 
by the drought treatment were up-regulated by the heat 
stress or recovery, such as VvbZIPS2 and 17 (Figure 7). 
This result suggests that different VvbZIP genes are in- 
volved in these two different biological pathways in re- 
sponse to drought and heat stress. 

Furthermore, the Pearsons correlation coefficient (PCC 
value) among pairs of VvbZIP genes based on all of the 
qRT-PCR expression data in the drought and heat stress 
conditions (Additional file 9) showed that 66 pairs of 
genes showed significant correlation at the 0.01 level 
(2-tailed), including 43 (65.2%) pairs showed signifi- 
cant positive correlation (PCC value ^0.99) and 23 



(34.8%) pairs showed significant negative correlation 
(PCC value g - 0.99). the fact that 17 positive and 13 
negative pairs of correlations between drought and 
heat stress conditionsuggests that some VvbZIP genes 
may also be involved in some common pathways in 
responses to drought and heat stresses. 

Discussion 

Increasing evidence suggests that the bZIP proteins play 
crucial roles in various developmental and physiological 
processes in plants, such as seed maturation [9,10], flower 
development [11-18], stress signaling [25,31,32] and patho- 
gen defense [30]. Much of the research into the functions 
of the bZIP genes have been characterized in model plants, 
especially in Arabidopsis, and some crops including rice, 
barley and Sorghum [2,4,7,54,55]. However, very little is 
known about this family in other species. 

In this report, we exhaustively identified the bZIP 
genes in grapevine genome, and examined their evolu- 
tionary history and high-quality expression profiles across 
the different developmental stages of various kinds of tis- 
sues and under drought and heat stress conditions. The 
information provides a comprehensive overview of the 
bZIP gene family in grapevine and lays the foundation for 
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future functional characterization of individual VvbZIP 
genes. 

The distribution and colinearity of grapevine bZIP gene 
family reflected evolutionary imprint of grapevine genome 

Grapevine is evolutionarily the earliest-diverging lineage of 
rosids and is a sister clade to eurosids I and II [36,42]. The 
species is believed to have evolved about 60 MYA, and its 
genome sequence revealed that it has not undergone any 
whole genome duplication since the gamma genome dupli- 
cation. Therefore, it serves as a good resource for under- 
standing the evolution and functional divergence of genes 
in rosids. The bZIP transcription factor family has evolved 
with a long history, which can be at least traced back to the 
earliest land plants, such as moss {Physcomitrella patens)^ 
which contains at least 37 bZIPs [55]. The grapevine, de- 
scribed as a palaeo-hexaploid; has homologous triplets in 
their genome and is believed to have resulted from a gen- 
ome fusion of a tetraploid and a diploid genome and subse- 
quent diploidization [35]; At the gene family level, the 
phylogenetic analysis of nucleotide binding site (NBS) gene 
family members [36], supported the fusion hypothesis, and 
suggested that the ancestral hexaploid genome of grapevine 
was derived from an allopolyploidization of two ancestral 
genomes: Va and Vc. However, the analysis of the calcium- 
dependent protein kinases (CDPK) gene family showed no 
indication of such imprint on fusion and divergence [56], 
and it may be because the gene family (17 members) is not 
large enough in tracing the evolutionary history, or the re- 
tention and loss after fusion was unequal. In this study with 
a large bZIP gene family (55 members), the distribution 
and gene duplication events of VvbZIPs (Figure 2), are very 
un-uniform. Five chromosomes (9, 10, 11, 16 and 17) har- 
bor no bZIP gene, and chromosome 15 harbored just one. 
Interestingly, chromosomes that lack bZIP genes seem to 
belong to the Vc genome, and this indicates that the grape- 
vine bZIP gene family suffered extensive gene loss after 
fusion and this loss mainly occurred in Vc genome. The 
12 homologous pairs of VvbZIP genes plus two triplets 
{VvbZIPO?/ 13/47 and VvbZIP14/22/37) occurred within 
the same homologous blocks and this partially supports 
the fusion hypothesis in the grapevine genome [35,36]. 

Furthermore, Comparing the number of VvbZIP genes 
with that of other species, it is -26% less than that in Ara- 
bidopsis, and -38% less than that in rice. This is probably 
because grapevine did not undergo any whole genome du- 
plication (WGD) events after the gamma WGD as shown 
in Arabidopsis, or the WGD after split between eudicots 
and monocots which happened in rice [57]. 

The colinearity-orthologues of the bZIP genes between 
grapevine and Arabidopsis 

Comparative genomic analyses across different taxa al- 
lows the transfer of functional information from a well- 



characterized taxon, such as the model plant Arabidop- 
sis, to another less-studied taxon, like grapevine. This is 
true with bZIP genes, many of which have been well- 
characterized in Arabidopsis, but few of which have been 
studied in grapevine. Therefore, the determining ortho- 
logous relationships between the Arabidopsis bZIP genes 
with those in grapevine (Table 1) will allow the grape- 
vine researchers to make a "best guess" of their functions. 
Most orthologues clustered at the same clade in the phylo- 
genetic tree (Figure 1) share common conserved motifs 
(data not shown). In Arabidopsis, AtbZIP36 (AREBl), 38 
{AREB2), and 37 {ABF3) have been shown to be master 
transcription factors that cooperatively regulate ABRE- 
dependent ABA signaling involved in drought stress re- 
sponse [25]. In grapevine, their orthologues are VvbZIPS 
and VvbZIP45, which were also shown to be highly indu- 
cible in response to drought stress in this study (Figure 7). 
Moreover, Nicolas et al. [39] recently demonstrated that 
VvbZIP45 transcripts mainly accumulated in the berry 
during ripening, and were up-regulated by ABA [39]. 
Comparative genomic analyses between grapevine and 
Arabidopsis provides a good reference for studying the 
biological functions of VvbZIPs. 

The grapevine bZIP gene family may participate in 
different aspect of organ development 

The in-depth analysis on a high-quality microarray data 
[50] (Figure 5) indicated that bZIP proteins may widely 
participated in the differentiation and development of 
many organs and tissues, and this was verified using qRT- 
PCR (Figure 6). Tissue-specific expression patterns include 
VvbZIP 19, 2S and S3 in seed development, VvbZIP29, 
32 and 54 in flower development, indicating their pri- 
mary roles in the morphogenesis and organ development. 
These results may lead to more directed understanding 
the function of these VvbZIP genes in grapevine develop- 
ment biology. 

Diverse VvbZIP genes are involved in responses to 
drought and heat stress 

Previous studies have shown that bZIP genes play an im- 
portant role in plants in response to various stresses, 
such as drought, salt, freezing temperature, and heat 
[2,4,46]. In this study, the expression profiles of VvbZIP 
genes under drought and heat stresses (Figure 7) re- 
vealed that grapevine bZIP genes are widely involved in 
responding to the drought and heat stress. Interestingly, 
some drought up-regulated genes were down- regulated 
in heat stress, which indicates that two sets of VvbZIP 
genes are involved in response to the drought and heat 
stress, respectively, others are co-expressed in response 
to both heat and drought stresses, and those genes may 
participate in shared roles between drought and heat 
stress. The data of the Pearsons correlation coefficient 
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(PCC value) among VvbZIP gene pairs also supports this 
notion (Additional file 9). 

Conclusions 

Grapevine {Vitis vinifera) genome contains 55 members 
in the bZIP gene family, and these genes are located in 
14 chromosomes. Genome distribution, gene organiza- 
tions and gene structures suggested a complex evolution- 
ary history of this family in grapevine. The grapevine bZIP 
family is greatly contributed by the segment/chromosomal 
duplications, which may be associated with the grapevine 
genome fusion events, their expression profiles in various 
tissue/organ and their development stages and responses 
to drought and heat stress condition demonstrated that 
this gene family is widely involved in grapevine organ de- 
velopment and drought- and heat-responses. The overall 
picture of this gene family and their potential involvement 
in growth, development and stress responses will facilitate 
further research on the bZIP gene family regarding their 
evolutionary history and biological functions. 

Methods 

Genome-wide identification of bZIP transcription factors 
in grapevine 

To identify all bZIP transcription factor genes in grapevine 
{Vitis vinifera) y all annotated proteins were downloaded 
from the latest database version (VI) of the 12X assembly 
of grapevine genome (http://genomes.cribi.unipd.it/). The 
Hidden Markov Model (HMM) profile of the bZIP domain 
(PF00170) downloaded from Protein family (Pfam) (http:// 
pfam.sanger.ac.uk/) was used for identification of the bZIP 
genes from the downloaded database of grapevine gen- 
ome using HMMER3.0 [58]. All output genes with 
default (< 1.0) E-value were collected and the online 
software SMART (http://smart.embl-heidelberg.de/smart/ 
set_mode.cgi?GENOMIC=l) was used to confirm the in- 
tegrity of the bZIP domain with E-value < 0.1 [59] and the 
incorrectly predicted genes were rejected. Finally, the 
non-redundant and confident genes were gathered and 
assigned as grapevine bZIP genes. Furthermore, the same 
methods were also used for detecting of possible add- 
itional bZIP signatures in other database, such as Geno- 
scope and NCBI, or identification of poplar bZIP genes. 
The protein sequences of grapevine, poplar and Arabidop- 
sis bZIP transcription factors were deposited at Additional 
file 10. 

Phylogenetic analysis and identification of conserved 
motifs of the VvbZIP family 

The conserved bZIP domain sequences of were extracted 
from bZIP proteins using online program SMART (http:// 
smart.embl-heidelberg.de/) for phylogenetic analysis. Mul- 
tiple alignments of protein sequences were executed by 
Clustal X 1.83 program [60]. Phylogenetic trees were 



constructed using MEGA 5.0 by the Neighbor- Joining 
(NJ), Minimal Evolution (ME) methods and the bootstrap 
test carried out with 1000 iterations, but only NJ tree was 
presented because all trees are consistent. The MEME 
program (version 4.8.1, http://meme.nbcr.net/meme/cgi- 
bin/meme.cgi) was used for identification of additional 
conserved motifs outside the bZIP domain, with the fol- 
lowing parameters: number of repetitions: any; maximum 
number of motifs: 50; and the optimum motif widths: 6- 
200 amino acid residues. 

Genomic structures, chromosomal locations and gene 
duplications for VvbZIP genes 

Each of the VvbZIP gene was positioned on grape chromo- 
some derived from the CRIBI (http://genomes.cribi.unipd.it/) 
using the programmed Perl script. The homeologous 
chromosome segments resulting from whole-genome 
duplication or fusion events were shown in the same 
color according to Jaillon et al. [35]. Gene structure 
display server program (GSDS, http://gsds.cbi.pku.edu. 
cn/index.php) was exploited to illustrate exon/intron 
organization for the bZIP genes by aligning the cDNAs 
with their corresponding genomic DNA sequences [61]. 
For detecting the gene duplication events, the MCScanX 
software (http://chibba.pgml.uga.edu/mcscan2/) was used 
and the E-value was set under lE-5. The other manipula- 
tions followed the operation manual. 

Plant materials preparation and treatment 

In vitro grapevine plants {V, vinifera, the sequenced geno- 
type PN40024) were maintained in vitro on MS 
medium supplied with 0.3 mg/L indole 3-butyric acid 
under a 16/8 h photoperiod (100 (imol m"^ s"^) at 25°C in 
the culture room. 

For drought stress treatment, one-month-old tissue 
cultured PN40024 plants were acclimatized in the 
lighted growth chambers [16-h light (25°C)/8-h dark 
(23°C), relative humidity of 60-70%, and light inten- 
sity of 20 (imol m-2 s-1 for 15 days, then watered regu- 
larly for 10 days, and this was set as the 0 d of treatment. 
One group of plants was then kept watered as the control, 
the other three groups were stopped being watered (the 
drought treatment) for a duration of 4, 8 and 12 days, 
respectively, and every duration of treatment had a cor- 
responding regularly-watered control. For each treatment, 
three independent plants (biological replications) were 
sampled. 

For heat stress treatment, in vitro-grown one-month old 
plants were directly placed in the 45°C lighted growth 
chamber (light intensity of 20 (imol m-2 s-1), excepting the 
temperature, the other conditions were the same as in the 
drought treatment. After heat treatment for 0, 2, 5 hours, 
plants were kept for 24 hours at 25°C for rapid recovery 
[27]. Each treatment contained three treated plants, or 
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biological replicates, and three corresponding non-treated 
control 

For tissue specificity expression of VvbZIP genes, tis- 
sues of root, stem, leaf (4-6* leaves from the tip) and 
shoot tips were harvested separately from three one- 
month-old tissue in vitro-grown plants of PN40024 
(three biological replicates). After freezing in liquid ni- 
trogen, all plant materials were stored at -80°C. 

RNA extraction, primer design and quantitative real-time 
PCR 

Total RNA was extracted using Huayueyang Quick RNA 
isolation Kit (Cat.No.: ZH120; Huayueyang biotechnol- 
ogy, Beijing, China), all protocols followed the manufac- 
turers procedure. To make sure that there have no 
DNA contamination, the traditional digestion protocol 
was used to eliminate DNA in the sample (rDNase®, Cat. 
No.: D2270A, TaKaRa Biotechnology, Dalian, China). The 
qualities and quantities of RNA were determined by agar- 
ose gel electrophoresis and Nanodrop ND-1000 Spectro- 
photometer (Thermo Fisher Scientific Inc.; USA). For 
cDNA synthesis, 500 ng high quality total RNA was re- 
verse transcribed with oligdT and random primers using 
SuperScriptlll Reverse Transcriptase (TaKaRa, Dalian, 
China) according to the manufacturer s instructions. 

Gene-specific primers for the VvbZIP genes were de- 
signed according to the predicted mRNA sequence using 
Beacon Designer 7.0 software (Premier Biosoft Inter- 
national, USA). To enhance specificity of the primer, 
most of them were cross-blasted the ORF and 3'UTR of 
VvbZIP genes, all of the primers were tested with PCR 
amplification, electrophoresis gels, melting curve ana- 
lysis and the products of PCR amplification were se- 
quenced at least three times, the maximum deviation of 
which were less than 3 base pairs. The sequences of 
primers and their products were listed in Additional file 
11, and noted that primers for VvbZIP28, 49, 52 and 53 
were not very specific, so it was just used in one or two 
treatment for determining the gene expression. 

For quantitative real-time PCR (qRT-PCR) assay, cDNA 
was diluted to about 100 ng/[A with ddH20. The qRT- 
PCR was carried out on ABI 7300 Real-time PCR System 
(Applied Biosystems) using SYBR® Premix Ex TaqTM 
(TaKaRa Code: DRR420A, TaKaRa, Dalian, China). The 
reaction system (total volume of 20 [A) was contained: 
10 \i\ SYBR" Premix Ex TaqTM (2x), 0.4 \i\ of each primer 
(10 (iM), 0.4 \i\ ROX Reference Dye (50x), 1 \i\ of template 
(about 100 ng/ (il), finally, add the ddH20 to total volume 
20 [iL The program were performed as the following steps: 
95°C/30 s for pre-denaturation (step 1), 95°C/5 s for 
denaturation (step 2), 60°C/34 s for primer anneaUng/ 
extension and gathering the fluorescent signal (step 3), 
then go to step 2 for 40 circles, each reaction was com- 
pleted in triplicate for the accuracy of results. At the end. 



the melting curve analysis was executed for verifying 
the specificity of the primer with the following program: 
95°C/15 s, 60°C/1 min, 95°C/15 s. Transcript levels were 
normalized against the average expression of VvACTIN 
gene (VIT_12s0178g00200), and the baseline and thresh- 
old cycles (Ct) were automatically determined by the soft- 
ware of the system. For the data analyses, 2'^^^ method 

[AACt=(Cttar get gene Ct Actin gene)treatment'' (Ct^arget gene 

CtActin gene)controi] was used for Calculating the relative ex- 
pression of VvbZIP genes. To visualizing the relative fold 
difference, all data were normalized based on setting up 
the relative expression level, the expression level of 0- 
point treatments for drought stress and heat stress and 
stem for plant tissue specificity was set as 1, "above 1" and 
"below 1" are considered as upper, or down, regulation in 
drought and heat treatment, and more, or less abundant 
in tissue specificity assay, respectively. The data (0 to 3 
folds) were displayed and clustered by MeV v4.8 software 
(http://www.tm4.org/), and the genes induced by each 
stress with their maximum folds (> 2 folds) of relative 
expression level were accounted and performed in the 
Figure 7. 

Microarray data analysis 

To understand the spatial and temporal expression patterns 
of VvbZIP genes during development, a high-throughput 
microarray data, from recent research [50], was employed 
for further analysis (Additional file 12). In the data sets, 54 
samples, including green and woody tissues and organs at 
different developmental stages as well as specialized tis- 
sues such as pollen and senescent leaves was determined 
(Additional file 6). The heatmaps were made by MeV v4.8 
software (http://www.tm4.org/) [62]. 

Statistics analysis 

To determine the co-expression relationship of VvbZIP 
genes in various tissues/organs development stages and 
drought and heat responses, the Pearsons correlation co- 
efficient (PCC) values among VvbZIP gene pairs were cal- 
culated based on the microarray data of various tissues/ 
organs development (Additional file 7) and qRT-PCR data 
of drought and heat stress (Additional file 9) using 
program Statistical Product and Service Solutions (SPSS 
v20.0). The significant correlation at the 0.01 level 
(2-tailed) was determined and denoted with red (sig- 
nificant positive correlation) and blue (significant negative 
correlation) background. 

Additional files 



Additional file 1: The gene structure of VvbZIPs. (a) The intron exon 
arrangement of VvbZIP genes; (b) the number of exons in VvbZIP genes. 

Additional flie 2: The position and pattern of introns within bZIP 
domain. 
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Additional file 3: Muitiieve! Consensus Sequence and their logo of 
VvbZIP proteins as predicted by MEME program. 

Additional file 4: The length of Leucine zipper region in bZIP 
domain of VvbZIP proteins. 

Additional file 5: Alignment of bZIP domain of VvbZIP proteins. 

Additional file 6: Description of organs and tissues of Vitis vinifera 
cultivar Corvina, from a vineyard in the Verona province (Montorio). 

Additional file 7: The Pearson correlation coefficient (PCC value) 
between VvbZIP genes based on the microarray data of various 
tissues/organs development. 

Additional file 8: Differentially expressed genes of grapevine bZIPs 
in various kinds of tissues. The microarray data was normalized based 
on setting up the expression level of first development stage of each 
tissue as 1, the folds of differential expression was subsequently calculated. 
The genes above and below the colored tissue bar indicated the 
up-regulated and down-regulated genes during this tissue development, 
respectively. The expression fold was presented at bracket along with each 
VvbZIP gene. The genes fluctuated more than 5 folds was highlight with 
red characters. 

Additional file 9: The Pearson correlation coefficient (PCC value) 
between VvbZIP genes based on the qRT-PCR data of drought and 
heat (45''C) treatment. 

Additional file 10: Primer sequences used for real-time PCR 
analysis. 

Additional file 11: Protein sequences of grapevine, poplar and 
Arabidopsis bZIP transcription factors. 

Additional file 12: High-throughput microarray data of VvbZIP gene 
expression. 
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